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ABSTRACT 

We present observations of the power spectral anisotropy in wave-vector space of 
solar wind turbulence, and study how it evolves in interplanetary space with increas- 
ing heliocentric distance. For this purpose we use magnetic field measurements made 
by the Helios-2 spacecraft at three positions between 0.29 and 0.9 AU. To derive the 
power spectral density (PSD) in , /c_L)-space based on single-satellite measurements 
is a challenging task not yet accomplished previously. Here we derive the spectrum 
PSD2D(k||, k_L) from the spatial correlation function CF2d(i'||) r^) by a transformation 
according to the projection-slice theorem. We find the so constructed PSDs to be dis- 
tributed in k-space mainly along a ridge that is more inclined toward the k_L than ky 
axis, a new result which probably indicates preferential cascading of turbulent energy 
along the k^ direction. Furthermore, this ridge of the distribution is found to gradually 
get closer to the k^ axis, as the outer scale length of the turbulence becomes larger 
while the solar wind flows further away from the Sun. In the vicinity of the ky axis, 
there appears a minor spectral component that probably corresponds to quasi-parallel 
Alfvenic fluctuations. Their relative contribution to the total spectral density tends to 
decrease with radial distance. These findings suggest that solar wind turbulence un- 
dergoes an anisotropic cascade transporting most of its magnetic energy towards larger 
k_L, and that the anisotropy in the inertial range is radially developing further at scales 
that are relatively far from the ever increasing outer scale. 

Subject headings: solar wind — turbulence — anisotropy 

1. Introduction 

Solar wind fluctuations are consi dered as the genuine and prominent example of magnetohy- 



drodynamic (MHD) turbulence (e.g.. iTu &: Marsch 



1995 



Goldstein et al. 1995: Bruno &: Car bone 
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20051 : iMarschI l2006l l . which is driven by solar activity and naturahy occurs in the inhomogeneous 
interplanetary space. Through in situ measurements made by the Helios 1 and Helios 2 spacecraft, 
the space plasma physics community has gained abundant knowledge about the radial evolution 
of solar wind turbu lence in the inner helio sphere. Magnetic power spectra of fluctuations at MHD 
scales were found (jBavassano et al.l Il982l ) to show two separate frequency ranges, whereby the 
power in the higher- frequency (f > 2.5 x 10^^ Hz in the spacecraft frame) range decays with radial 
distance (like r~^'^) faster than that in the lower- frequency (f < 2.5 x 10~'^) range (with radial 
scaling like r~'^'^). 

The radial evolution of the lower-freq uency inagnetic power spec tra can be reproduced by the 
WKB-theory of Alfven wave propagation (|Whangj|l973l : lHollwegill974l ). which predicts a similar ra- 
dial evolution. Whereas the higher-frequency magnetic power spectra, which show a steeper profile 
(Kolmogrov-like) with its spectral break frequency shifting toward s lower values during th e radial 
evolution, were successfully reproduced by Tu's turbulence model (jTu et al.lll984l : iTulllQSa ). which 
took into account (together with the WKB description) the nonlinear interaction between counter- 
propagating imbalanced Alfven waves. Moreover, the n ormalized cross-helic i ty (Alfvenicity) was 



shown to decrease w ith increasing heliocentric distance (jRoberts et al 



1987 



Marsch fc Tu 



199C 



Grappin et al.lll990l ). which to explain was beyond the scope of Tu's model. To self-consistently 
describe the radial evolution of turbulent energy, cross-helicity, and Alfven ratio, su bstantial theo- 
retical efforts had to be made, which finally resulte d in general transport equations ([Marsch fc Tu 



198S 



Tu Marsch 



1990 



Zhou Sz Matthaeuslll989l ) for the related spectra. 



In numerical simulations of MHD turbulence, the assumed background magnetic field (Bq 
was found to cause spatial anisotropy of the turbulent fluctuations along and across the mean field 
with the parallel scale generally being larg er than the perpendicular scale (e.g 



Shebalin et al 



1983 



Biskamp &: Miilleri l2000l : ICho et al.l |2002| ) . For balanced strong MHD turbulence with vanishing 

2 /3 

cross-helicity, the anisotropy is predic ted to reveal a scaling rela tion obeying ky ~ k^ , which was 
derived in a phenomenological theory ( Goldreich fc Sridharil 19951 ) based on the conjecture of critical 
balance, i.e. the rough equality between the linear wave-propagation time and nonlinear eddy- 
interaction time. Numerical simulations further showed that balanced strong turbulence behaves for 
strong or weak Bq differently in its scaling properties across Bp : Iroshnikov-Kraic hnan scaling was 
found for strong Bq and Goldreich-Shridar scaling for weak Bq (iMiiller et al.ll2003l ). These different 
scalings are argued to be probably attributed to an increase of dynamic alignment as the cascade 
pro ceeds to smaller scale, which may also induce scaling anisotropy in the plane perpendicular to 
Bo (iBoldvrevlboOsI ). 



However, in the solar wind and particular in fast streams one usually observes imbalanced 
turbulence with outgoing waves dominating over incoming waves. This imbalanced turbulence 
implies different nonlinear interaction time scales for the oppositely propagating waves, and is 
more complex than the balanced one . Its physical nature remains a controversial issu e , although 



several theories have been proposed (iLithwick et al.l 120071 : iBeresnyak LazarianI l2008l : IChandran 
20081 : IPodesta Bhattacharieelboid ). 



- 3 - 



The spatial anisotropy of solar wind tur bulence was studied by means of data analysis employ- 



ing various tools, such as c orrelation function (iMatthaeus et al 



2009l : IWicks et alJl201(3l . 



2011 



20071 : iNemeth et alJl20ld). structure function (|Luo k Wull2010l : IChen et al.ll201Q. i2012l) power of 



1990 



magnetic components (IBieber et alj 119961 ) . and power law scaling (jHorburv et alj 120081 : 



Dasso et a l. 2005: jOsman &: Horburv 



Podesta 



). Scaling anisotropy becomes more clearly visi ble if one uses a scale - 
dependent local mean magnetic field (Bojocai) (which was first introduced bv lHorburv et alJ (120081) 



apply ing the wavelet technique) instead of a constant global mean magnetic field (jXessein et al 



20091 ). Some efforts have been also made to recon struct magnetic PSD i n multi-dimension al wave- 
vector space by means of the k- filtering method ([Sahraoui et al.l l20ld : iNarita et al.l |2010| ) , which 
was develo ped originally to distingu ish a limited number of plane waves from multi-position mea- 



surements (Pincon &: Lefeuvre 



1991 



y). The integrated PSDiD(k_|_) with a spectral index ~ —1.6 as 
obtained from k-filtering method seems to resemble the reduced PSDiD(k_|_) as derived from the 
direct wavelet transformation. However, in some c ases studied with the k-filtering method, the inte- 
grated PSDiD(k||) shows a spectral index ~ —7.0 (jSahraoui et al.ll2010l ) ^ differing significantl y from 



the reduced PSDi£)(k||) (~ k,, ) as obtained by the wavelet method ( Horburv et al.l l2008l ) . The 



reliability of k-filtering method for estimating turbulent power spectra may need further validation, 
e.g. by applying it to numerically simulated turbulence with known scaling. 

Previous studies have revealed the evolution of reduced ID-PSD in the inner heliosphere, and 
have presented evidence of wave- vector anisotropy at specific positions [e.g., 1 AU]. However, the 
turbulence anisotropy pattern at 0.3 AU (the innermost distance reached in-situ so far) and its 
evolution trend between 0.3 and 1.0 AU has not yet been investigated. To do this is an important 
task, because it will provide the needed information about the evolution of the energy cascading 
route in k-space, and reveal possible ways of turbulent energy dissipation required for sustained 
solar wind heating. This work is dedicated to a study of MHD turbulence anisotropy and will 
provide new knowledge on its spectral characteristics. The data analysis to achieve these goals is 
briefly described as follows. 

Firstly, we estimate the second-order structure function as a function of 0rb (the angle between 
the radial direction and the local mean magnetic field vector). Accordingly, the angular distribution 
of the spatial correlation function is obtained, using the relation between structure function and 
correlation function. Secondly, we fit the measured structure function with a compound fit function, 
resembling a power-law dependence at short scale and giving an exponential trend at large scale. 
The fitted angular correlation function is subsequently derived. Thirdly, under the assumption of 
a statistically time-stationary state, the relative 2D-PSD in (k||,k_|_) space is constructed from the 
fitted angular distribution of the correlation function, whereby we make use of the pro j ection -slice 
theorem which is fundamental for image processing in medical tomography (see iBovikI (|2000l ) for a 
detailed review). 
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2. Analysis method 

In this section, we describe the applied methods and the data analysis, which includes: how 
to derive angular distributions of the structure function SF(r, ^vb) and the correlation function 
CF(t, 6'vb); how to fit SF(r, ^vb) and CF(r, ^vb) appropriately; and how to obtain PSD2D(k||, k_|_) 
as transformed from CF(r, ^vb)) which is in turn obtained from CF(r||,r_L) by assuming a quasi- 
steady state with ry ~ KswTcos^vb and r± ~ KwTsin^vB, corresponding to the Taylor assumption 
of fluctuations being frozen into the flow, and thus being simply convected by the wind past the 
spacecraft. Here V"sw is the solar wind speed. 

The second-order magnetic structure function is defined as the ensemble average of the squared 
magnetic field vector difference. It can be written as 

SF(r) = ((B(t + ^)-B(t-^))2), (1) 

where the angular bracket denotes in practice a time average in our subsequent data analysis. 
This time average permits one to quantify the global scaling of the magnetic fluctuations, without 
distinguishing a possible scaling-law difference for different angles (^vb) between the sampling 
direction and the local mean magnetic field vector (Bq, local)- The local mean magnetic field is known 
to be changing in time and depend on scale (Bo,iocai(^) ''"))) leading to a scale-dependent variation 
of the angle ^VB with time. For solar wind with a radial speed much larger than the velocity 
fluctuation amplitude, the quantity ^VB can be approximated by 0rb (i-e., the angle between the 
radial direction and the Bojocai direction), which is used hereafter. To estimate the structure 
function value at a certain time scale r' and for a certain angle 9'-^^, one needs to pick out the 
values of SF(i,r') at those times when ^RB(i)''"') = ^rb> and then make an average over all the so 
picked samples. Therefore, the corresponding angular distribution of the structure function can be 
expressed as 

, , /o^(B(t + ^)-B(t-4))^^t|,,,,,,)^,,^ 
Jo dt eRB{t,r')=eJiB 

Here the time period for the whole chosen data set is indicated as T. It should be much larger than 
the time scale r', and thus we may formally take the limit T oo. 

Expressing the ensemble average used in equation ([1]) explicitly as a time average, the relation 
between the structure function SF(t) and the correlation function CF(t) can be obtained from the 
subsequent calculation: 



£ B\t + ^)dt + £ B\t - ^)dt - 2 £ Bit + ^) ■ Bit 



1 
T 

2CF(r = 0) - 2CF(r). 



(3) 
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Using the above definition ([2]), the angular distribution of the correlation function can also be 
approximated by the angular distribution of the structure function, yielding on the basis of ([3]) the 
following relation: 

SF(t, e') = -2 (B{t + I) . B(t - ^)) \e,,=e' +{B^it + ^)) \e^,=e' + {b^I - ^)) 

~-2CF(r,e') + 2CF(T = O,0') (4) 
~ -2CF(r,e') + 2CF(t = 0), 

where angular isotropy of CF at r = was assumed in the derivation. Under Taylor's hypothesis 
the solar wind fluctuations can be considered time stationary, as the wave phase speed is small in 
comparison to the supersonic convection speed, and then CF(t,9') can be rewritten as a spatial 
correlation function in the 2D r-space, 

CF(r,0~CF2D(r||,r^), (5) 

with ry = VswTCOsO' and r± = Vsv^t s'mO' . This completes the derivation of the two-dimensional 
correlation function from the structure function. We note that the frozen-in-flow Taylor's hypothesis 
may be slightly weakened for smaller heliocentric distance with smaller Alfven Mach number, which 
drops from higher than 10 at 1 AU to 3-4 near 0.29 AU. The quantity of main interest is the 
power spectral density PSD2D(k|| ; kx), which in can in principle be obtained directly from Fourier 
transformation of CF2d(i"|| ) r^) as follows: 

/ + 00 /■ + 00 
/ CF2D(r||,r^)exp(-i(k||r||+kxr^))dr||dr^. (6) 
-oo J — oo 

However, we take here a new route to estimate PSD2D(k|| , k^). It can also be derived fro m the 
proje cted (integrated) ID correlation function on the basis of the projection-slice theorem (jBovik 



2000) with help of the following formula: 



/ + 00 /• + 00 
/ CF2D(r||,rx)exp(-i(k(r||Cos0k + r^sin0k)))dr||drx 
-oo J — oo 

/ + 00 /• + 00 
/ CF2D(r'cos^k-u'sin0k,r'sin0k + u'cos0k)exp(-i(kr'))dr'du' (7) 
-oo J — oo 

■oo 

CFiD(r';0k)exp(-i(kr'))dr', 



where 6k is the angle between k and Bo,iocah and CFiD(r';0k) is the ID projection (integration) of 
CF2D(r||;r±) along the direction normal to k, 

/ + CO 
CF2D (r' cos 0k - u' sin ^k, r' sin 0k - u' cos 0k)du'. (8) 
-oo 



Therefore, there are two approaches to calculate PSD2D(k|| , k^), one may adopt either Equa- 
tion [6] or [71 In practice, the estimation of the 2D correlation function with help of Equation [6] 
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introduces some uncertainty, as the noise involved in the data may destroy the required positiv- 
ity of the PSD in the entire {k\\,k±) space. To guarantee this positivity of PSD everywhere, one 
needs to approximate the CF with some kind of positive-definite fit function before the Fourier 
transformation. It is hard to find an adequate function that globally fits the observed CF2d(i'|| , r_L) 
well, whereas it is relatively easy to choose a proper fitting function for the projected CFi£)(r; 0^). 
Therefore, in our work, we will use a fitted CFiD(r;6'k) to reconstruct reliably the PSD2D(k||, k±) 
according to Equation [Tl 

To provide the reader with an intuitive impression about the relations between CF2d(i'||, r_L), 
CFiD(r, ^k), and PSD2D(k|| , k_|_), we present the schematic illustration shown in the upper panel 
of Figure [H which explains the two roads from CF2D to PSD2D (direct 2D Fourier transform and 
indirect method based on the projection-slice theorem). Similarly, one slice of CF2D at certain 
angle 6r is also the ID inverse Fourier transform of PSDid as projected from PSD2D onto the 
corresponding direction k with = 9r, an approach which is displayed in the lower panel of 



Figure [H The relation betwee n CF^n and PSE> 



was used in previous studies (jMatthaeus et al 



n is t he basic method for calculating CF^n, whic h 



1990; iDasso et al.ll2nn,'5l : lOsman fc Horburvl I2nn7l ) . 



In principle, it is also possible to derive PSD2D from PSDid according to the method of inverse 
Radon transform (filtered back-projection) (private communication with M. Forman). However, 
this method fails in a typical benchmark test due to extreme large PSD at small |k|, which blurs 
the entire reconstructed PSD2D thereby destroying its original pattern. 

Speaking of the fitting function for the CF, we need to mention also the fitting function for the 
SF, which is used to reproduce the key features of the SF. For example, people usually adopt an 
exponential function to fit the profile of the SF at large scale, while they use a power-law function 
for the small-scale trend. However, as far as we know, there exists no attempt to describe both 
the small-scale power-law trend and the large-scale exponential trend simultaneously with a single 
fitting function. To fulfill this task, we suggest a compound function, 

SF(r) = 2Ro-[l-exp(-(-)P)], (9) 

which interpolates between these limits. There are three parameters to be fitted: Rq means the 
auto-covariance at r = 0, represents the correlation time at large scale, and the index p describes 
the power-law scaling at short scale. Generally, for SF(t, 0rb) at different ^rb, the parameters Rq 
and Tc do not change a lot, while p remains variable. Therefore, in our practice, i?o and Tc are 
obtained by fitting the time-averaged SF(r), and then p is determined at various 0rb by fitting 
SF(t, 0rb), but only after Rq and Tc were set. Another practical reason for presetting Rq and Tc 
before fitting SF(t, 0rb) is that for every 0rb the calculated SF(r, 6'rb) is usually unable to reach 
to the outer scale. 
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3. Data analysis results 

The magnetic data (with a time resolution of about 0.25 s) used here is from measurements by 
Helios-2 spacecraft at three radial positions (0.29, 0.65, and 0.87 AU) during three time intervals 
(day of year: 106-109, 76-78, and 49-51 in 1976). The solar wind streams explored during these 
time intervals are known to be recurrent streams emanating from a common source region on the 



Sun (iBavassano et al.lll982l ). The corresponding radial evolution of ID reduced magnetic PSD was 



presented in that paper , which observa tionally promoted the development of the WKB-like solar 



wind turbulence model (|Tu et al.lll984l ). Three decades later, we analyse the same data set again, 
but for the purpose of revealing the evolution of solar wind turbulence in terms of its wave- vector 
anisotropy. 

We use Equation [2] to estimate the second-order structure function SF(r). It is defined as 
the magnetic vector difference squared, which is averaged respectively over the three time intervals 
of our data set. During the estimation, the data gaps are excluded without making any type of 
interpolation. The top three panels of Figure [2] illustrate the estimation results as red curves. The 
blue lines are fitting results based on Equation [U which basically match the estimates at both 
small and large scales. The fitting parameters (i?o[nT^], tc[s], p) at three radial positions are found 
to be: (827, 116, 0.61), (53, 465, 0.61), and (25, 857, 0.67), respectively The fitting parameter 
tc (corresponding to the correlation time) increases with heliographic distance. The values of the 
exponent p relate to the power-law index (~ — (p + 1)) of the corresponding PSDs, which is found 
to be around —1.6, i.e. near the Kolmogorov value of —5/3. The bottom three panels of Figure [2] 
show the corresponding correlation function CF(r) as derived from Equation [3l 

We calculate the structure functions in the angular dimension as a function of 0rb according 
to Equation [2l and display them in the first row of Figure [3l Apparently, the distribution of SF is 
not uniform in the angle range between 0° and 90°, with a lower level near 0°. The non- uniform 
angular distribution is more significant at short scales [e.g., < 100 s]. For SF(r, ^rb) at larger scales 
(r > 100 s), it gradually changes from uniformity at 0.29 AU to non-uniformity at 0.87 AU. This 
angular non-uniformity is a feature hinting at anisotropy of the power spectrum in the wave- vector 
space. Likewise, the extension of the angular non-uniformity towards larger scales indicates that 
the wave-vector anisotropy of larger-scale fluctuations evolves as heliocentric distance increases. 

We also fit the estimated structure function by the function SF(r, ^rb) of Equation [9l To make 
sure the fitting process converges for every angle, we restrict the number of fitting parameters to 
p, while we fix the other two parameters (i?o and tc), both of which may be regarded as constant 
without angular dependence. The fitted angular distributions are illustrated in the second row of 
Figure O which look similar to the observations. The angular dependence of the fit parameter p 
is plotted in the third row, showing that the angular variation of SF(r, ^rb) is non-uniform not 
only in magnitude (first row in Figure [3]) but also in the scaling index (third row) . We note that 
SF(t, 0rb) as shown in Figure [3] relates to the squared module of the magnetic- vector difference 
{SB"^ + 5By + 6B1). The structure function SF(r, 0rb) for the component 6B?, (parallel to Bo,iocai) 
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shows a similar non-uniform angular dependence. However, the calculated SF for (5i3^ has a plain 
segment starting at small r, and cannot be fitted well by the function of Equation [9l 

The angular distribution of the correlation function CF(r, 0rb) is derived from the fit function 
SF(r, 0rb) according to Equation [H In the light of the projection-slice theorem as applied to 
the relationship between the 2D functions CF and PSD (lower panel in Figure [1]), the quantity 
CF(r, ^rb) is essentially a 2D correlation function CF(r|j , r_|_), which is in principle an inverse Fourier 
transform of the 2D PSD(k|j,k_L) yet not known. In Figure HI we plot the resulting CF(r||,r_L). 
The coordinates of the abscissa (ry) and ordinate {r±) are estimated by ry = 1/sw''"cos0rb and 
= Kiw'T" sin 0RB ) respectively. The main part of CFfrii , r i ) is elongated alo ng rn, which is similar 



to the "2D" population of the so-called Maltese cross (jMatthaeus et al.lll99Cll ). However, the "slab" 



population, which was reported in previous stat istical studies of CF with r ii parallel to th e direction 



of interval-averaged (non-local) magnetic field (jMatthaeus et al.lll990l : iDasso et al.ll2005l ). is not so 
prominent in our cases. 

Ideally, the corresponding PSD(k|| , k_|_) can be gained directly from 2D Fourier transform of 
CF2D(r||! r±)- However, in practice, the transformed value might be negative or not certainly posi- 
tive, thereby restraining the application of the direct 2D Fourier transform. To obtain PSD2D(k|| , k±), 
we then turn to Equation [7] for a step-by-step derivation. Firstly, by integrating CF2D(r||) ^x) over 
the path normal to the direction with certain angle 6' with respect to ry, the reduced ID CFiD(r) 
corresponding to the angle 9' is calculated. Secondly, the corresponding PSD as a Fourier trans- 
form of CFiD is calculated. To guarantee the positivity of the estimated PSD, CFid is fitted before 
transformation with a function related to that for SF as previously described. According to the 
projection-slice theorem, the estimated PSD profile is essentially a slice of PSD2D(k||, k_|_) along k 
with 6' with respect to Thirdly, PSD2D(k||, k_|_) is formed by assembling various PSD profiles, 
with different angles ranging from 0° to 90° with respect to fcy. We note that, in calculation, CF2D 
and CFiD one cannot let r go to infinity. As a result, the transformed PSD2D may slightly depart 
from the real one. Therefore, in Figure O we just present the normalized PSD2D,n rather than the 
absolute PSD2D- The uncertainty (confidence interval) for the estimated PSD2D is not provided 
here, since due to the complexity of the estimation method that was not yet possible. 

Obviously, the normalized PSD2D,n shown in Figure[5]is not uniformly distributed at all angles, 
indicating an anisotropic wave-vector distribution. This anisotropy is mainly characterized by a 
ridge distribution which has a bias towards k± as compared to k^^. Moreover, as the heliographic 
distance increases, the ridge distribution becomes more inclined toward k± at the same |fc|, in 
association with lower PSD (darker blue in the figure) around the /cy region and higher PSD 
(brighter blue in the figure) around the k± region. The discovery of this bent ridge and its radial 
evolution imply that solar wind turbulent energy cascades preferentially along the k± as compared 
to the /cy axis, and the turbulence cascade radially develops with more energy cascading to the k± 
region, as the scale (1/|A;|) is shifting away from the radially-growing outer scale (l/|A;o|). In addition 
to the major ridge distribution, a minor population seems to exist close to fcy (see Figure Eh.), and 
appears to become weaker at farther distances (see Figures [5)3, c). The observational fact that PSD 
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is composed of two populations, with the major one bending more perpendicularly and the minor 
one becoming weaker, seems compatible with the previously suggested two-component turbulence 
model, which invokes non-damping convective structures (spatiall y varying across B p) that are 
superposed on damping Alfven waves (spatially varying along Bq) (|Tu Marschlll993l ). 



To emphasize the trend of the ridge distribution and its radial evolution, we estimate the 
ridge position of every scale by averaging the angles with local Ig(PSD) as the weights (i.e., first- 
order moment centroid method). The estimated ridge positions are shown as black dashed lines 
in Figure [5] and appear straight. Whether or not the straightness is realistic is yet unknown. 
Furthermore, we fit the estimated ridge position with following simple formula. 



1/3 ,2/3 



(10) 



where A;o(= 27r/(Kjw "^c)) is related to the outer-scale wave-number, a is the coefficient to be fitted, 
which is ~ 3.2, 3.9, and 3.9 for our three cases. We find that the simple /cy — k_|_-relation profile 
according to Equation [10] is basically coincident with the observed ridge distribution. However, some 
departures, e.g., the estimated black dashed line looks more straight than the fitted red line, still 



remain. Nevertheless, the relation 



ky^k^/^)^ as predicted by the critical-balance hypothesis 
(jGoldreich Sridhaiill995l ) for MHD turbulence, seems to describe well the observed anisotropy of 
solar wind turbulence. The role of A:o, which was once neglected in previous observational studies, 
in shaping the anisotropy shall be emphasized here. It may be the reduction in which causes 
the development of the spectral anisotropy (increasing inclination toward at the same \k\) in 
interplanetary space as heliographic distance increases. 

Solar wind heating mechanism may be inferred from the radial evolution of the ridge trend. 
According to linear Vlasov theory, Alfven waves with plasma /3p € [0.1, 1.0] usually become dis- 
sipated due to proton cyc lotron resonance when t hey have fcyc/wp G [0.1,1.0], where c/uj^ is 
the proton inertial length (|Garv Nishimural 12004 ) . On the other hand. La ndau resonance be- 
comes more and ni ore prominent as plasma j5p rises (iGary &: Nishimurall2004l ) and k^pg increases 
(jHowes et al.l 120061 ) . where pg is the proton gyroradius. iHowea (j201lh pointed out that, Landau 
damping calculated in the gy ro-kinetic limit is not sufficient for the empirically estimated proton 
heating (jCranmer et al.ll2009l ) at small heliocentric distances {R < 0.8 AU). At these small dis- 
tances, turbulent cascade is spec ulated to ap proach to proton cyclotron frequency before being 
terminated by Landau resonance (lHowesll201ll ). However, the relative contributions from cyclotron 
resonance and Landau resonance to solar wind heating at different radial distances have not yet 



been addressed from observations. The approximated relation (/cy = a 



k 



1/3 ,2/3 







k/ with a G [3,4]) 



for the observed ridge distribution may be used to address this issue. In Figure [U we just simply 
extend the approximated ridge profile in a larger wave- vector space to see what kind of resonance 
(cyclotron or Landau) would probably terminate the cascade. As a result, we find that, at the 
three distances within 1 AU, the extended ridge profile clearly approaches to cyclotron resonance 
(marked by k^^c/up > 0.5) ahead of Landau resonance (marked by kj_pg > 1.0). Moreover, one may 
expect that, as the distance increases (> 1 AU) in association with reduction of k^, the ridge profile 
would first exceed the threshold k±pg = 1 before approaching to k^<^c/up = 0.5, which implies a 
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dominance of Landau damping over cyclotr on damping at larger dis tances. According to critical- 
balance theory in MHD and kinetic regimes (jScliekochihin et al.ll2009l ). the extension of ridge in the 
MHD inertial range may be still related to Equation \TU[ while the extension part in the kinetic 
(dissipation) range may deviate from Equation [10] with more inclination towards k± . For sake of 
simplicity, we neglect such deviation of extended ridge in the kinetic range from that in the inertial 
range. 

The value of coefficient a is also worth emphasizing here. If a were one third of the ap- 
proximated value (3.3/3 = 1.1), the ridge profile (red dash-dot-dot line in Figure [6|) would exceed 
k±pg = 1 without approaching to k^^c/ujp = 0.5, leading to insuffi cient h e ating rate by Landau res- 
onance within 1 AU according to the gyro- kinetic prescription by iHowed (|201ll ). On the other side, 
if a were too large (saying 3.3 x 3 = 10), the ridge profile (red dashed line in Figured]) would lie well 
below k±pg = 1, implying the absence of transiti on frorn cyclot ron damping to Landau damping 
around 1 AU (inconsistent with the conclusion by ]Howes] ()201l] )l. Therefore, the a value (€ [3,4] 
obtained here) besides kQ is another important parameter for grasping the essence of solar wind 
heating mechanism, a may be expressed as the ratio of e to V^ko with e being the energy cascade 
rate if fc y = (e/Vl)^/^k'^/\ which is usually assumed in critical-balance theory ()Goldreich &: Sridhai] 



199,4 ]Schekochihin et al.l]2009l ^ 



4. Summary and discussion 

We have made the first successful attempt to reconstruct, on the basis of single spacecraft 
measurements, the 2D spectral density PSD2D(k||) k^) for solar wind MHD turbulence. We esti- 
mate the angular distribution of the second-order structure function SF(r, ^rb)) and derive the 
corresponding correlation function CF2D(r|j > r±)) which in principle is an inverse 2D Fourier trans- 
form of PSD2D(k|| ) k^). The transformation from time scale r to spatial scale r, when building 
up CF2d(i'||; r_|_), is based on Taylor's hypothesis that solar wind fiuctuations are quasi-stationary 
within the flow transit time scale, as the solar wind passes by the spacecraft. The 2D direct Fourier 
transform of CF2D(r||, rx) fails to guarantee the required positivity of PSD2D(k|| , k^). Alterna- 
tively, we employ for the first time a method based on the projection-slice theorem, which connects 
the integrated CFid{t,6') with the corresponding slice PSD2D(k, ^') of the PSD via a ID Fourier 
transform, to fulfill that task. Before the ID Fourier transformation, CFiD(r,^') is fitted smoothly 
to guarantee the positivity of the transformed PSD2D(k, 0')- 

As a result, SF(r, ^rb) shows a non-uniform angular distribution with more power being located 
in the perpendicular region (^rb ~ 90°) than in the parallel region (0rb ~ 0°) of wave- vector space. 
Moreover, there is angular dependence of the scaling law for SF(r, ^rb) at short scales, whereby the 
scaling index p drops from ~ 0.9 at ^b = 0° to ~ 0.6 at 0rb = 90°. We find that SF(r, Orb) have 
at all three positions (0.29, 0.65, and 0.87 AU) the two above properties, indicating the prevalence 
of anisotropy in the turbulence throughout the inner heliosphere. This result obtained within 1 AU 



-Il- 



ls similar to that found for the SF anisotropy beyond 1 AU (|Luo &: Wull20ld ). The corresponding 
correlation functions CF2D(r||,r_|_) clearly show that magnetic fluctuations are correlated at longer 
(shorter) length along (across) the background magnetic field. 

The corresponding PSD2D(k|| , k^) at the positions within 1 AU is revealed to have a ridge dis- 
tribution with a bias towards k± as compared to fcy, suggesting a preferential cascading along k±. 
This kind of ridge distribution has never be en reported in prev i ous studies at 1 AU, e .g. those based 
on the wave-telescope (k-filtering) method (jNarita et al.ll201(]| : ISahraoui et al.ll2010l ). Furthermore, 
this ridge distribution is found to become ever more inclined toward the k± axis with increasing 
heliographic distance, thus indicating a radial development of the wave-vector anisotropy. The ob- 
served radial evolution of the ridge casts new light on the scaling relation between /cy and k±, which 

1/3 2/3 

may empirically be approximated by fcy ~ ak^ kj^ , with a € [3, 4] and k^ being the wave-number 
of the outer scale. This approximation for the wa ve-vector anisotropy seems to indicate critical- 
balance- type cascading ([Goldreich &: Sridhaiill995l ) of solar wind turbulence. A possible influence 
of kQ on the anisotropy development, which was neglected in previous observational analyses, is 
also found. 

However, the evolution of the ridge distribution cannot represent the whole story about wave- 
vector anisotropy of solar wind turbulence. There seems to be a minor population located near 
which is beyond the scope of critical-balance turbulence theory. The apparent two-component 
distribution of PSD2r)(k||, k i ) seems to be conne cted with previous two-component models, e.g., 
models wit h "slab" -|- "2D" iMatthaeus et al.l ll99d). models composed of Alfven waves and convected 
structur es (ITu &: MarschI 119931). and conject ures with critical-balanced component plus slab com- 



ponent (jForman et al 



2011 



He et al.ll2012bl ). We find that the minor population seems to weaken 



further with increasing heliographic distance, leaving more energy distributed in the region close 
to the k± axis. This gradual migration of energy towards k± might indicate a relative enrichment 
of turbulence energy carried by convective structures and explain the observed assoc iated shortage 
of Alf venicity, which was already discussed in the previous two-component model by lTu fc: Marsch 
(| 19931 ). The observed "sl ab" -like minor coinponent is crucia, l for scattering o f energetic particles in 



the interplanetary space (jBieber et al.lll996l : 



Chandranll200d : iQin et al.ll2002l ). The radial evolution 



of anisotropic turbulence may be quantified in the future and incorporated into the transport model 
of energetic particles. 

The estimated PSD2D(k|| , k^) is believed to im pose valuable observational constrai nts on the 
theoretical models of solar wind turbulence. Recently. I Cranmer &: van BallegooijenI (120121 ) modelled 
PSD2D(k|| , k^) at different heliocentric distances by solving a set of 2D cascade-advection-diffusion 
equations, with the total power pre-determined by the damped wave-action conservation equa- 
tion and the reduced PSDiD(k^) pre-set by the ID advection-diffusion equation. Their modelled 
PSD2D(k|| , k^) looks partly similar to our observational spectrum, in the sense of where the major 
power is located. However, the differences in distribution pattern and radial evolution between 
observational and modelled spectra call for a substantial improvement of the models for solar wind 
turbulence. 
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The approximated ridge profile as extended to large k\\ and may give a hint about the 
resonance type responsible for solar wind heating at different radial distances. The extended ridge 
profile at small distances {R < 0.8 AU) is found to reach larger fcy where proton cyclotron reso- 
nance acts before Landau damping sets in. As the distance increases, the extended ridge profile, 
which is inclined more towards k±_ due to the reduction of feg, tends to arrive at Landau resonance 
before cyclotron resonance. Such performance of the approximated ridge profile confirms obser- 
vationally previous conjecture about t he transition from cyclotron resonance to Landau resonance 



with increasing heliographic distance (jHowesI 12011 



Our results are just limited to the MHD inertial range, but the analysis should be extended 
to kinetic scales where several typical properties have been revealed: steeper power-law mag- 



netic spectrum (Sahraoui et al 



Salem et al 



2011 



2012; 



200S 



Alexandrova et al. 



(IBale et al.l 12005 ). enhanced mag netic compressibility (jSmith et al.l l2006l : [Hamilton et al 



Podesta &: Gary 



2011 



20091'). enha.nced electric-field spectrurn 



He et al.ll2012al). and two-c omponent pattern in the magnetic helicity (jHe et al 



He et al. 



2012 



30). 



200S: 



These observations seem to be in favour of the 
oblique Alfven waves or kinetic Alfven waves (KAW) as the candidate for explaining the dominant 
fluctuatio ns in ion-scale turbulence. T he oblique Alfven/ion-cyclotron waves may be via resonance 



diffusion (jMarsch &: Bourouaindl201lh responsible for the formation of the observed wide proton 



beam. The theory of KAW itsel f and its role in kinetic turbulence have been studied intensively 



(IHoUwedllQQgl: 



Howes et al 



Wu fc Chao 



2011 



2004 



Howes et al. 



20081 : 



Zhao et al 



2011 



Voitenko &: de Kevser 



2011 



). Ion cyclotron waves, which are con sidered responsible for the ion p erpendicular 



heating (IBourouaine et al.ll2010l ). were also identified (jjian et ahlbood : Ihc et~alll201ll ). There are 
other possible wave modes, e.g. fast wh istler waves , ion Bernstein waves, and fast-cyclotron waves , 
which may exist in kinetic turbulence (iGary et al.l l2012l : iTenBarge et al.l 120121 : IXiong &: Lil l2012l ) . 
Spectral break at the ion-kinetic scale seems to be almost constant (about 0.5 Hz in the space- 
craft frame) with radial distance (jPerri et al.ll2010l : IBourouaine et al.ll2012l ). The spectral break 
frequency might corresponds to the proton inertial length in quasi-2D turbulence when considering 
a large-scale backgrou nd magnetic field Bp (w hich is obtained through averaging over a time period 
higher than 1 hour) ([Bourouaine et al.ll2012l ). However, not much is presently known about the 
radial evolution of solar wind turbulence at ion-kinetic scales. 

In the future, with the help of high-time-resolution measurements to be made by the wave 
and particle instruments flown on such mission like Solar Orbiter and Solar Probe Plus, the radial 
evolution of the wave-vector anisotropy at kinetic scale may be studied, and more new results will 
be obtained on the spectrum anisotropy in the inertial range that was analysed here. 
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Fig. 1. — Sketch of relation between 2D-PSD and 2D-CF based on the projection-shce theorem. 
(Top) Two approaches to derive PSD2D(k|[,k_|_) from CF2d(i'||i rx): direct 2D Fourier transform 
and indirect Fourier transform of the projected CFiD(r, 6). (Bottom) Vice versa for the derivation 
of CF2D(r||,rx) from PSD2D(k|| , k^) 
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Fig. 2. — (Top) Time averaged second-order structure functions based on Equation [T] (red) at three 
positions (0.29, 0.65, and 0.87 AU) and their corresponding fitting results according to Equation [9] 
(blue). The sets of the three fitting parameters {Rq [nT^], tc [s], and p) are (827, 116, 0.61), (53, 
465, 0.61), and (25, 857, 0.67) at 0.29, 0.65, and 0.87 AU. (Bottom) Corresponding correlation 
functions based on Equation [3] (estimations in red, fitting results in blue). 
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Fig. 3. — (Top) Angular distribution of second-order structure functions SF(r, 0rb) estimated on 
the basis of Equation [2j (Middle) Fitting results for SF(r, 0rb) with SF(t) at every 0rb being 
fitted according to Equation [H whereby p is fitted in angular dependence. (Bottom) Fit parameter 
p as a function of 0rb, revealing the scaling anisotropy of the structure function. The error-bars 
denote the fitting errors of the parameter p. 
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Fig. 4. — Angular distribution of CF2D as derived from SF(t, 0rb) according to Equation U] and 
displayed in (r||,r_|_) space under the Taylor hypothesis of near time-stationarity. An elongation 
of CF along ry, implying the location of most turbulent energy close to k±, is revealed at all 
positions. The proton gyroradius pg is 17, 48, and 70 km at 0.29, 0.65, and 0.87 AU. It is used for 
normalization of the spatial coordinates. 
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Fig. 5. — PSD2D(k||, k_|_) at three positions (0.29, 0.65, and 0.87 AU) as derived from CF2d(i'|| , r^) 
according to Equation [7] following the projection-slice theorem. The major components (ridge 
distribution with its centroid position aligned as black dashed line) may be roughly described by 
Equation [To] shown as red solid line, which means the wave- vector anisotropy develops as the outer- 
scale wave-number (ko) becomes smaller with increasing heliocentric distance. A weakening trend 
of the minor component that is inclined to ku becomes also visible. 
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Fig. 6. — Implication of solar wind heating mechanism from the extension of ridge profile. Red 
solid lines denote the extended ridge profiles = a ■ ■ k^/^ with a = 3.2, 3.9, 3.9 for the three 
panels) as superposed on PSD2D(k|| , k±) in larger wave-vector space. Large A;|| with k^^c/cop > 0.5 
(to the right of blue dashed line) indicates the region of ion cyclotron resonance. Landau resonance 
becomes active when k±pg > 1 (above magenta dashed line) and dominant for k^^c/up <^ 1.0. Radial 
evolution of the approximated ridge profile and its intersections with the threshold lines {k^^c/uip = 
0.5 and k±pg = 1) indicate the transition of cascade termination from cyclotron resonance to Landau 
resonance as the solar wind flows further away. Ridge profiles with larger (smaller) a{= 10(1.1)) (red 
dashed (red dash-dot-dot) lines) would lead to cyclotron resonance (Landau resonance) separately 
at all distances without a transition, which seems unrealistic. The proton gyroradius (proton 
inertial length) pg{c/ujp) is 17(42), 48(82), and 70(121) km at 0.29, 0.65, and 0.87 AU. 



